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Abstract 

Background: We investigated the relationship between average monthly temperature and the most common clinical 
pathogens causing infections in Intensive care patients. 

Methods: prospective unit-based study In 73 German intensive care units located in 41 different hospitals and 31 different 
cities with total 188,949 pathogen isolates (102,377 Gram-positives and 86,572 Gram-negatives) from 2001 to 2012. We 
estimated the relationship between the number of clinical pathogens per month and the average temperature in the month 
of Isolation and in the month prior to isolation while adjusting for confounders and long-term trends using time series 
analysis. Adjusted incidence rate ratios for temperature parameters were estimated based on generalized estimating 
equation models which account for clustering effects. 

Results: The incidence density of Gram-negative pathogens was 15% (IRR 1.15, 95%CI 1.10-1.21) higher at temperatures 
>20°C than at temperatures below 5°C. E. cloacae occurred 43% {IRR=1.43; 95%CI 1.31-1.56) more frequently at high 
temperatures, A. baumannil 37% (IRR = 1.37; 95%CI 1.11-1.69), S. maltophilla 32% (IRR =1.32; 95%CI 1.12-1.57), K. 
pneumoniae 26% (IRR = 1.26; 95%CI 1.13-1.39), CItrobacter spp. 19% (IRR =1.19; 95%CI 0.99-1.44) and coagulase-negatlve 
staphylococci 13% (IRR = 1.13; 95%CI 1.04-1.22). By contrast, S. pneumoniae 35% (IRR = 0.65; 95%CI 0.50-0.84) less 
frequently isolated at high temperatures. For each 5''C Increase, we observed a 3% (IRR= 1.03; 95%CI 1.02-1.04) increase of 
Gram-negative pathogens. This increase was highest for A. baumannil with 8% (IRR= 1.08; 95%CI 1.05-1.12) followed by K. 
pneumoniae, CItrobacter spp. and E. cloacae with 7%. 

Conc/usion: C\\n\ca\ pathogens vary by incidence density with temperature. Significant higher Incidence densities of Gram- 
negative pathogens were observed during summer whereas S. pneumoniae peaked In winter. There is increasing evidence 
that different seasonality due to physiologic changes underlies host susceptibility to different bacterial pathogens. Even If 
the underlying mechanisms are not yet clear, the temperature-dependent seasonality of pathogens has Implications for 
infection control and study design. 
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Introduction 

Seasonality is characteristic of many infectious diseases and has 
been recognized since Hippocrates who began his first book 
Epidemics with the description of the weather [1]. He wrote in 400 
BC: "Whoever wishes to investigate medicine properly should 
proceed thus: in the first place to consider the seasons of the 
year...". 

A strong seasonal effect can be seen in many respiratory or 
bacterial gastrointestinal infections and in seasonally recurring 
childhood infections [2,3]. Respiratory infections follow seasonal 
patterns. In temperate climates, respiratory illness is most common 
in winter months, and in tropical settings the incidence of lower 



respiratory-tract infections is generally increased during the rainy 
season [4]. 

By contrast, Acinetobacter infections show infection rates twice as 
high in late summer months as in early winter months [5] . Data on 
summer peaks in the incidences of other Gram-negative bacterial 
infections among hospitalized patients have only recently been 
published [6,7,8]. The authors of the study by the University of 
Maryland reported that incidence rates for E. coli were 1 2 percent 
higher in summer than in winter. Rates for E. cloacae were even 46 
percent higher with the same seasonal variation. Seasonal 
variation has been described also for non-infectious diseases like 
depression, heart attack or stroke [9,10]. 
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Multicenter data on the seasonality of Gram-positive and 
Gram-negative clinical isolates in ICUs remain scare [11,12]. 
Thus, the aim of our study was to look for temperature 

associations of pathogens in a network of geographically variant 
intensive care units in Germany and to discuss reasons for and 
implications of seasonal variation. 

Materials and Methods 

Ethics Statement 

All data were anonymous and collected in accordance with the 
German recommendation for good epidemiological practice with 
respect to data protection [13]. As a federal law, the German 
Protection against Infection Act (Infektionsschutzgesetz §23) 
regulates the prevention and management of infectious disease 
in humans. All hospitals are obliged to continuously collect and 
analyse continuously healthcare-associated infections and drug- 
resistant pathogens [14]. These data are reported regularly to the 
National Reference Centre for the Surveillance for Nosocomial 
Infections. Ethical approval and informed consent were therefore 
not required. 

Intensive Care Units 

From January 2001 to December 2012, a total of 73 German 
adult intensive care units (ICU) reported data to SARI (Surveil- 
lance of Antimicrobial Use and Antimicrobial Resistance in 
Intensive Care Units) [15]. The ICUs are located in 41 different 
hospitals and 34 different cities in Germany. 

Data for air-conditioning were available for 64 of the 73 ICUs 
(88%). Thereof, 83% were air-conditioned during the entire 
analysis period, 6"/o were air-conditioned during parts of the 
analysis period and 1 1 % were not air-conditioned. 

Data Collection 

Monthly susceptibility data were collected from microbiology 
laboratories for 5 Gram-positive pathogens {S. aureus, Coagulase 
negative staphybcocci (CoNS), E.faecalis and faecium, S. pneumoniae) and 
8 Gram-negative pathogens {E. coli, P. aeruffnosa, K. pneumoniae, E. 
cloacae, S. maltophilia, S. marcescens, Citrobacter spp., A. baumannii). Copy 
strains - defined as an isolate of the same species showing the same 
susceptibility pattern throughout the period of one month in the 
same patient, no matter what the site of isolation - were excluded. 
During data collection, it was not possible to differentiate culture 
data from isolates obtained less than or equal to 48 hours after 
ICU or hospital admission from isolates obtained greater than 48 
hours (nosocomial). Monthly patient and device days were 
collected from ICU personal. Temperature data were collected 
as means per month in cadi individual year for cities or states from 
the publicly available database of the Deutscher Wetterdienst [16]. 

Statistical Analysis 

The incidence rate of isolates (IR) was defined as being the 
number of isolates of a pathogen per 1000 patient days. We 
investigated the following temperature parameters: 

1 . temperature in the month of isolation (categorized by 5 degree 
Celsius (°C) steps) to investigate the association with temper- 
ature including the proportionality (linearity) of the association; 

2. temperature in the month prior to isolation (categorized by 5°C 
steps) to investigate the temporal shift/ delay in the association 
(to reflect endogenous colonisation), as well as the proportion- 
ality of the association. 



Incidence rates with 95% confidence intervals were calculated 
for each pathogen for the 2 above-mentioned parameters. 

Time series analysis using generalized linear models were 

performed to estimate the association of the number of pathogens 
per month with temperature parameters. 

Since obserx^ations within an ICU are not statistically indepen- 
dent due to the diagnostic policies (especially the frequency of 
microbiological tests), adjusted incidence rate ratios (IRR) with 
95% confidence intervals (CI) were estimated. They were based on 
generalized estimating equation (GEE) models which account for 
this clustering effect by using an exchangeable correlation 
structure [17]. 

We used negative binomial distribution in the models instead of 
Poisson distribution because the variance exceeds the mean and 
because overdispersion was observed for all type of pathogens. The 
Lagrange multiplier test was used to test whether the negative 
binomial model significandy differs from the Poisson model. The 
log number of patient days during each month was treated as an 
offset in the model. 

The following parameters considered confounding in the 
frequency of the occurrence of pathogens: long-term trends 
between 2001 and 2012 (linear, quadratic and cubic time trend), 
the number of pathogens in the previous month as an 
autoregressive factor (ARl), and the following structural param- 
eters: type of ICU (medical, surgical, interdisciplinary), size of ICU 
(</> = 12 beds), type of hospital in which ICU is located 
(university, teaching, other) and size of hospital in which ICU is 
located (</> = 1000 beds). AH confounding parameters were 
parameterized as continuous or dummy parameters and added 1 
degree of freedom to the model. 

The Multivariable Model Building Strategy was 
Performed in Two Steps for each Pathogen 

1 . All confounding parameters were included in a full GEE 
model, and non-significant parameters were excluded stepwise. 
The selection criterion was the smallest Chi-square value and 
p> = 0.05 in the type III score statistic. 

2. For each pathogen, a GEE model was calculated with the 
significant confounding parameters from the first step and the 
temperature parameters. Temperature parameters were con- 
sidered the categorical variable in the model with four degrees 
of freedom in order to test differences in the frequency of 
pathogens between temperature intervals; temperature inter- 
vals were considered the ordinal/linear covariate in the model 
with one degree of freedom in order to test the linearity of the 
association. 

An additional analysis was performed with the following 
monthly process parameters: length of stay (LOS, da)"s); and 
device use (device days per 100 patient days) for central venous 
catheters, urinary tract catheters and ventilators to adjust for the 
severity of illness. This additional analysis showed the same 
significant results as without these parameters. However, the 
process parameters were only available in 85 percent of months 
analyzed. Therefore, we present our results without adjusting for 
LOS and device use. 

The quasi-hkelihood information criterion (QIC) as a modifi- 
cation of the Akaike information criterion (AC I) was used as 
goodness-of-fit measure in the GEE model. P-values less than 0.05 
were considered significant. All analyses were performed using 
SPSS [IBM SPSS statistics, Somer, NY, USA] and SAS [SAS 
Institute, Gary, NC, USA]. 
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Table 1. Baseline characteristics of 73 German ICUs, January 2001 to December 2012. 





Parameters 


Statistic 


Structure parameters of ICU and hospital 


Number ICUs {number hospitals; number cities) 


73 (41; 34) 


Type of ICU 


interdisciplinary, No. (%} 


34 (46.6) 


medical, No. (%) 


18 (24.7) 


surgical. No. (%) 


21 (28.8) 


Type of hospital in which ICU is located 


university. No. (%) 


23 (31.5) 


teaching, No. (%) 


40 (54.8) 


other. No. (%) 


10 (13.7) 


Size of ICU (No. beds), median (IQR) 


12 (10-18) 


Size of hospital in which ICU is located {No. beds), median (IQR) 


854 (488-1229) 


Analyzed months 


months total. No. 


5965 


Months per ICU", Median (IQR) 


84 (36-127) 


Distribution of temperature C'C) 


Monthly temperature, median (IQR); (Min; Max) 


9.9 (4.4-16.2); (-5.6; 24.4) 


Process parameters" 


Incidence rate of pathogens'^ {No. pathogens/1000 patient days), median {IQR) 


78 (61-97) 


Length of stay and device use (N =62 ICUs) 


Length of stay {days), median (IQR) 


4.0 (3.2-5.5) 


Ventilation utilization (No. ventilation days/100 patient days), median (IQR) 


43 (31-57) 


Central venous catheter utilization {No. central venous catheter -days/100 patient days), median {IQR) 


74 (57-84) 


Urinary tract catheter utilization (No. urinary tract catheter-days/1 00 patient days), median {IQR) 


85 (76-91) 



Abbreviations: No., number; ICU, intensive care unit; IQR, interquartile range; ''C, Degree Celsius, 
"calculated per ICU. 

"^considered 5 Gram-positive pathogens (S. aureus, Coogulase negative staphylococci, E. faecalis and faecium, S. pneumoniae) and 8 Gram-negative pathogens {£ coli, P. 
aeruginosa, K. pneumoniae, E. cloacae, S. maltophilia, S. marcescens, Citrobacter spp., A. baumannii}; Copy strains were excluded. 
doi:l 0.1 371 /journal.pone.0091 1 0S.tOOl 



Results 

73 ICUs reported data on 2,243,856 patient days (pd) and 
188,949 pathogen isolates (102,377 Gram-positives and 86,572 
Gram-negatives). Table 1 shows the characteristics of the 73 ICUs. 

In median, 84 months per ICU and 78 pathogens per 1000 
patient days were analysed fromjanuary 2001 to December 2012. 
The median monthly temperature was 9.9°C and varied between 
-5.6 and 24.4°C. 

The time series of the incidence rates for the different Gram- 
positive and Gram-negative pathogens in 73 German ICUs and 
mean monthly temperature in Germany are depicted in figure 1 
and figure 2. Generally, the incidence of the Gram-negatives 
follows the increase in the temperature and follows a decrease in 
temperature likewise. This is not the case for Gram-positives like. 
S. pneumoniae or S. aureus. 

Table 2 gives the number and the incidence rate of pathogens 
by temperature in the month of isolation over the 12-year period, 
stratified by 5°C temperature interval steps. 

Temperature in the Month of Isolation 

The adjusted incidence rate ratios (IRR) of the temperature in 
the month of isolation are shown in figure 3 and in more detail in 
table 3. 



In total, the incidence density of Gram-negatives pathogens was 
15% (IRR 1.15, 95%CI 1.10-1.21) higher in months with high 
temperatures (^20°C) when compared to months with low 
temperatures (<5°C). 

In eight of the thirteen pathogens (two of the five Gram- 
positives; six of the eight Gram-negatives), the adjusted IRR for 
temperature in the current month was significant (type III test). In 
S. pneumoniae, the IRR was 0.65 (95%CI 0.50-0.84) for months > 
20°C compared to months <5°C. In other words, S. pneumoniae 
occurred 35 percent less frequently at high temperatures (&20°G) 
than at low temperatures (<5°C). 

Figure 1 shows that the incidence density of S. pneumoniae and 
the temperature were in opposite phases: when it was cold, the 
incidence of S. pneumoniae was high. 

The other seven pathogens were significandy positively associ- 
ated with high temperatures (&20°C). E. cloacae occurred 43% 
(IRR= 1.43; 95%CI 1.31-1.56) more frequendy at high temper- 
atures, A. baumannii 37% (IRR =1.37; 95%CI 1.11-1.69), S. 
maltophilia 32% (IRR= 1.32; 95%CI 1.12-1.57), K. pneumoniae 26% 
(IRR= 1.26; 95%CI 1.13-1.39), Citrobacter spp. 19% (IRR= 1.19; 
95%CI 0.99-1.44) and coagulase-iiegative staphylococci 13% 
(IRR= 1.13; 95%CI 1.04-1.22). 
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Figure 1. Time series of montKily aggregated incidence densities 
blue line) in 73 German ICUs, January 2001 to December 2012. 

doi:10.1371/journal.pone.0091105.g001 

Linearity 

Testing the linearity of the relationship (testing temperature 
parameters as ordinal/linear covariate in the model with 1 degree 
of freedom) shows that the relationship was linear (figure 3, table 3). 
For the Gram-negative pathogens, the IRR was 1.03 (95% CI 
1.02-1.04) per 5°C step. For each 5°C increase, we observed a 3% 
increase in the incidence density of Gram-negative pathogens. 
Specifically, there was an 8% increase for A. haumannii 
(IRR=1.08; 95%CI 1.05-1.12), foUowed by K. pneumoniae, 
Citobacter spp. and E. cloacae with 7%. There was a 9% decrease 
for S. pneumoniae (IRR = 0.91; 95%CI 95%CI 0.88-0.94). 

Temperature in the Month Prior to Isolation 

To account for the temporal shift between the temperature and 
the frequency of pathogens, we considered the temperature in the 
month prior to isolation. Table 4 gives the number and the 
incidence density of pathogens over the 12-year period, stratified 
by 5°C steps temperature intervals in the month prior to isolation. 

The adjusted incidence rate ratios (IRR) for the temperature in 
the month prior to isolation are shown in figure 4 and in more 
detail in table 5. 

In eleven of the thirteen pathogens, the adjusted IRR for 
temperature in the month prior to isolation was significant. For S. 
aureus, E. faecium and S. pneumoniae, the association was negative 
compared to month prior to isolation <5°C; for the seven Gram- 
negatives (excluding S. marscecens) and CoNS, it was positive. 
Moreover, testing the linear associations showed significance 



Coagulase negative Staphylococci 
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Gram-positive patKiogens (bold red line) and temperature (thin 



(excluding E. faecalis and 6*. marscecens). Moreover, the differences 
were more pronounced in comparison to the results which take the 
temperature in the month of isolation into account. 

Discussion 

The four main findings of this study in 73 ICUs and on 188,949 
non-duplicate bacterial cultures over a 12-year period are as 
follows: 

i. Eleven of the thirteen analysed pathogens were associated 
with temperature iS. pneumoniae, E. coli, K. pneumoniae, E. 
cloacae, S. maltophilia, Citrobacter spp., A. baumannii, CoNS; S. 
aureus, E. faecium and P. aeruginosa). Only E. faecalis and S. 
marcescens showed no temperature association. 

ii. The Gram-positive pathogen CoNS and the Gram-negatives 
E. coli, P. aeruginosa, K. pneumoniae, E. cloacae, S. maltophilia, 
Citrobacter spp., A. baumannii were positively associated with 
temperature. In other words, the warmer the weather, the 
higher the incidence density. S. aureus, E. faecium and S. 
pneumoniae were negatively associated with temperature; i.e. 
the colder the weather, the higher the incidence density. 

iii. The associations were proportional to the temperature 
intervals (or inversely proportional, respectively). 

iv. The association of the frequency of isolation and tempera- 
ture for Gram-negatives was stronger when the mean 
monthly temperature in the month prior to isolation was 
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Figure 2. Time series of monthly aggregated incidence density of Gram-negative patKiogens (bold red line) and temperature (thin 
blue line) in 73 German ICUs, January 2001 to December 2012. 

doi:10.1371/journal.pone.0091105.g002 



considered instead of tlie temperature in the month of 
isolation. 

Explanations for these differences are mainly hypothetical. 
There is conflicting evidence in support of different ideas, 
suggesting that there is no one single reason but rather a collection 
of different factors. One reason could be that environmental 
factors like temperature, light, humidity and wind influence the 
frequency of pathogen isolation. However, room temperature and 
other environmental factors are usually fairly constant in ICUs, in 
hospitals and in laboratories. Therefore, it seems more likely that a 
variation in the endogenous flora of the patient plays a role. 

Most infections, including nosocomial infections, are endoge- 
nous in nature. In a comprehensive study in German ICUs, only 
14.5 percent of nosocomial infections were associated with 



exogenous patient-to-patient transmission [18]. The mostly 
endogenous origin of infection or colonisation could also explain 
why the association of frequency of occurrence with temperature a 
month before the isolation of the pathogen was the same or even 
more prominent than the association with the current outside 
temperature. This holds true for infections that are generally 
classified as community acquired, such as infections caused by S. 
pneumoniae, as well as those mostly classified as nosocomial like 
infections by A. baumannii. But irrespective of the endogenous 
origin, why do infections with some Gram-positive and Gram- 
negative bacteria peak in different seasons of the year? 

The phenomenon that pneumococcal infections increase in 
winter has been repeatedly documented [19,20]. The incidence 
density ratio of S. pneumoniae was 35 percent lower in month with 
temperatures > = 20°C (compared to months <5°C) in our study 
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Figure 3. Adjusted incidence rate ratios for Gram-positive and Gram-negative pathogens by temperature in the month of isolation. 

As reference we defined temperatures in the month of isolation <5"C. For the calculation we used generalized estimating equations (GEE) models 
using negative binomial distribution and logarithmic transformed patient days as offset parameter. Data based on 5965 monthly data sets of 73 
different German ICUs, January 2001 to December 2012. The adjusted models accounts for cluster effects within an ICU and for significant 
confounding parameters from the first model building step: S.aureus coagulase negative staphylococci E. faecalis £ faecium ^' S. 



pneumoniae ^\ E. coli P. aeruginosa K. pneumoniae '"° <:°"f°""der)^ ^ cloacae 



, 5. maltophilia ^' S. marcescens ^, Citrobacter 



spp. '"° ^°"<°""'''^'\ /\, baumannii * in which: ^ autoregressive factor AR(1) number pathogens in the previous month; ^ linear time trend; 

^ quadratic time trend; ^, cubic time trend; ^ ICU is located in a university hospital; hospital size > = 1000 beds; ^, type of ICU medical; ^ type of ICU 
surgical. *P<0.05 in chi-square test with 4 degrees of freedom (type III test) in the adjusted GEE model; Whiskers represent 95% confidence intervals. 
doi:1 0.1 371 /journal.pone.0091 1 05.g003 



as well, which is just the opposite of what we found in Gram- 
negatives. Dowell et al. evaluated data of seven geographic areas in 
the United States and found that pneumococcal disease correlated 
inversely with temperature [19]. One might expect higher rates 
throughout the year if colder temperature is what drives the winter 
increase. Paradoxically, the coldest states had the lowest rates in 
their study. Dowell believes that the photoperiod (dark-light cycles) 
and not temperature per se contributes to seasonality. The 
incidence of infections may therefore be due to changes in the 
susceptibility of the host to the particular pathogen. Susceptibility 
might differ for different pathogens. It might depend on 
photoperiod being the most pervasive signal for seasonal changes 
in the biological system and may be mediated by the length of the 
daily melatonin pulse (melatonin, a pineal hormone, has been 
shown to mediate seasonal adjustments in immune function) and a 
range of other physiologic parameters like the serum Cortisol level, 
IL-6, CD4 and CD8 ceUs or circulating B ceUs [21]. 

Other explanations could be pathogen appearance and 
disappearance or host behaviour changes like seasonal variation 
in the consumption of food or outdoor activities [21]. The periodic 
occurrence of gastrointestinal infections and outbreaks in summer 
is well recognized [22]. However, the objection can be made that 
the most prominent summer peak of Gram-negative clinical 
isolates was by A. baumannii, a pathogen which does not appear to 
be a typical environmental organism but has its main reservoir in 
humans. By contrast, other Acinetobacter species indeed seem to be 
distributed widely in nature. A. calcoaceticus, for example, is found in 
water and soil and on vegetables [23]. 

In a landmark study, Perencevich et. al. analyzed 26,624 unique 
chnical cultures over 8 years in a US-American university hospital 
and observed significant summer peaks in the incidences of non- 
duplicate clinical cultures of E. coli, A. baumannii, E. cloacae and P. 
aeruginosa [6] . This was true for Gram-negative but not for Gram- 
positive pathogens. Interestingly, the results of our study were very 
similar to Perencevitch et al, who found a 46 percent increase of E. 



cloacae in summer compared with winter. In our study, there was 
43 percent increase in months with temperatures S20°C 
(compared to months <5°C). There was a 21 percent increase 
for A. baumannii and 37 percent respectively and a 12 percent 
increase for E. coli and 5 percent respectively. Perencevitch et al 
could not show a seasonal variation for S. aureus, but found that 
Enterococcus species were 9 percent less frequent in summer than in 
winter. The authors concluded that the different findings for 
Gram-positive and Gram-negative pathogens reduced the likeli- 
hood for confounding variables such as hospital-related factors, 
and underline the assumption that there are mechanistically 
different reasons for summer peaking only by Gram-negatives. In 
an other study, 211,697 blood stream infections were analysed in 
132 U.S. hospitals over 8 years. They looked at the monthly 
outdoor temperature and showed that they were associated with 
substantially increased frequency of bloodstream infections among 
important Gram-negative bacteria [8]. These results are in line 
with our findings: an increase in monthly temperature of 5.6°C 
corresponded to independent increase in A. baumannii by 10.8% 
versus 8% our study, E. coli 3.5 versus 1%, K. pneumonia 8 versus 
5% and P. aeruginosa 7.5 versus 3%. However, they did not 
investigate Serratia, Stenotrophomonas, Enterobacter and Citrobacter spp. 

Our analysis considered hospital-related confounding parame- 
ters as well (long-term trends, number of pathogens in the previous 
month, and structural parameters, among others). 

Other studies did not see seasonal variation in bloodstream 
infections (BSI) caused by Enterobacter or Klebsiella. These studies 
were limited by a very low number of pathogens and focus on one 
single pathogen, however [24,25]. Anderson et al. investigated 
whether K. pneumoniae, Enterobacter or Serratia species had higher 
incidences during the warmest month of the year in four hospitals 
on four continents [26]. They found that increasing environmental 
temperature was independently associated with an increased 
incidence of K. pneumoniae BSI but not for Serratia or Enterobacter 
BSI. The authors explained their findings by the fact that K. 
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Figure 4. Adjusted incidence rate ratios for Gram-positive and Gram-negative pathogens by temperature in the month prior to 
isolation. As reference we defined temperatures in the month prior to isolation <5 C. For the calculation we used generalized estimating equations 
(GEE) models using negative binomial distribution and logarithmic transformed patient days as offset parameter. Data based on 5965 monthly data 
sets of 73 different German ICUs, January 2001 to December 2012. The adjusted models accounts for cluster effects within an ICU and for significant 
confounding parameters from the first model building step: S.aureus ''' coagulase negative staphylococci 
pneumoniae £ coli P. aeruginosa K. pneumoniae '"° ""'"-^"der)^ ^ cloacae ^' ^' S. maltophilia 

spp. '"° ™"'°""der)^ ^ baumannii * in which: ^ autoregressive factor AR(1) number pathogens in the previous month; ^ linear time trend; 

^ quadratic time trend; ^, cubic time trend; ^ ICU is located in a university hospital; hospital size > = 1000 beds; ^, type of ICU medical; ^ type of ICU 
surgical. *P<0.05 in chi-square test with 4 degrees of freedom (type III test) in the adjusted GEE model; Whiskers represent 95% confidence intervals. 
doi:10.1371/journal.pone.0091105.g004 
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b 



pneumoniae is the most heat tolerant of all enteric pathogens and 
that Klebsiella survives better at higher humidity. 

In contrast with Anderson et ai, our study revealed a linear and 
prominent temperature dependency oiE. cloacae, We also found no 
seasonality for S. marcescens, which is in accordance with their 
findings. It remains unclear why S. marcescens seems to differ from 
other Gram-negatives. S. marcescens can also survive and grow 
under extreme conditions but the epidemiology is somewhat 
different from other Enterobactereaceae: Serratia appears to be less 
likely to colonize the gastrointestinal tract and more likely to 
colonize the respiratory or urinary tracts of hospitalized patients 
[27] [28]. The review by Leekha et al. on the seasonality of S. 
aureus showed that some studies could show an increase in 
incidence with increasing temperature but not an increase in 
colonizations. Leekha et al. determined, however, that there are 
important gaps in the literature on the seasonality of S. aureus 
infections, the most evident of which being a lack of appropriate 
statistical methodology (i.e. use of time-series analysis or time series 
regression models) for seasonahty assessment, or the studies' being 
performed in a single center and/ or a single year [29] . By contrast, 
our study showed a small but significant decrease of S. aureus 
isolates of 1 % (infections and colonizations) with a temperature 
increase of 5°C. In contrast to other pathogens with seasonal 
effects, this is of little clinical relevance. 

Generally, an improved understanding of the forces driving 
seasonal physiologic changes and endogenous biologic rhythms 
may uncover unimagined opportunities for disease prevention. 
However, it was beyond the scope of our study to investigate 
physiology and susceptibility of Gram-positive and Gram-negative 
pathogens in different seasons. Regardless of the mechanism 
responsible for infection, our multicenter data have contributed to 
the emerging evidence for the seasonality of Gram-negative 
bacteria, and might have an impact for infection control, therapy 
and future study designs. Estimating the magnitude of the impact 
of those parameters provides more detailed information with the 



aim to prevent infections. The knowledge that Gram-negative 
infections were on average 15 percent more likely to occur in 
summer compared with winter could help in the selection of 
optimal empirical antibiotic therapy. Studies should consider 
seasonality as additional parameter. For example, studies - even 
randomized controlled trials - analyzing infections caused by 
Gram-negatives in summer versus in winter could considerably 
bias the results if they do not control for seasonality. 

The strengths of our study are the multicenter design with 73 
ICUs in 41 hospitals and 34 different cities in Germany, the 
inclusion of 13 Gram-positive and Gram-negative pathogens 
(n= 188,949 non-duphcate bacterial cultures), the relatively long 
study period of 12 years and the consideration of all isolates, i.e. 
not only isolates from blood culture. Additionally, our time-series 
regression models took cluster effects within an ICU into account, 
and considered different confounding parameters and adjusted for 
severity of illness of patients (length of stay and device use within 
the ICU). 

Our study has also several limitations. First, all isolates were 
included and we could not differentiate by localization, origin 
(nosocomial or not regarding hospital and/or ICU), colonization 
or infection. Second, conclusions from our study may not be 
generalized for other settings because we analyzed only adult ICU 
patients. Third, our epidemiological data were obtained from an 
ongoing surveillance system which is ward-based, not patient- 
based, therefore making patient-specific analyses impossible. 
Forth, we could not exclude pathogens from outbreaks or 
screening which might have inflated the number of pathogens 
isolated. 

It is interesting to speculate that seasonal variations in 
physiologic parameters might underlie the distinct seasonal 
variation of Gram-positive or Gram-negative pathogens by 
altering the susceptibility of host cells to infection. Even if the 
underlying mechanisms are not yet clear, the temperature- 
dependent seasonality of Gram-positive or Gram-negative path- 
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ogens has implications for infection control (surveillance, choice of 
empirical antibiotic treatment) and study design. 
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